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Abstract 

A review is presented on all higher order QCD corrections to deep inelastic structure 
functions. The implications of these corrections for polarized and unpolarized deep inelastic 
lepton-hadron scattering will be discussed. 



To appear in the proceedings of the 1996 HERA Physics Workshop 



1 Introduction 



The past twenty years have shown much progress in the field of perturbative calculations in 
strong interaction physics [1]. This in particular holds for the radiative corrections to the deep 
inelastic structure functions. Sometimes these corrections could be even extended up to third 
order in the strong coupling constant a s . The structure functions we would like to discuss are 
measured in deep inelastic lepton-hadron scattering 

h{k 1 ) + H{p)^l 2 {k 2 )+"X" (1) 

where h, l 2 stand for the in- and outgoing leptons respectively. The hadron is denoted by H and 
"X" stands for any inclusive hadronic state. The relevant kinematical and scaling variables are 
defined by 

q = k 1 -k 2 q 2 = -Q 2 >0 x = $- y = 2L (2) 

2pq pki 

with the boundaries 

< y < 1 < x < 1 (3) 

Reaction (1) proceeds via the exchange of one of the intermediate vector bosons V of the standard 
model which are represented by V = 7, Z, W. In the case of unpolarized scattering with V = 7 
one can measure the structure functions Fl(x,Q 2 ) (longitudinal) and Fi(x,Q 2 ) (transverse) or 
the better known F 2 (x, Q 2 ) which is related to the former two via 

F 2 (x, Q 2 ) = 2x F l (x, Q 2 ) + F L (x, Q 2 ) (4) 

When V = W or V = Z one can in addition to Fi, F 2 and Fl also measure the structure function 
F 3 (x,Q 2 ) which is due to parity violation of the weak interactions. In the case the incoming 
lepton and hadron are polarized one measures besides the structure functions Fj (i = 1,2, 3, L) 
also the spin structure functions denoted by gi(x,Q 2 ) (i — 1, • ■ -5). At this moment, because 
of the low Q 2 available, reaction (1) is only dominated by the photon (V = 7) so that one has 
data for gi(x, Q 2 ) (longitudinal spin) and g 2 (x,Q 2 ) (transverse spin) only. The measurement of 
the structure functions at large Q 2 gives us insight in the structure of the hadrons. According 
to the theory of quantum chromodynamics (QCD) the hadrons consist out of quarks and gluons 
where the latter are carriers of the strong force. When Q 2 gets large one can probe the light cone 
behaviour of the strong interactions which can be described by perturbation theory because the 
running coupling constant denoted by a s {Q 2 ) is small. In particular perturbative QCD predicts 
the (Revolution of the deep inelastic structure functions mentioned above. Unfortunately the 
theory is not at that stage that it enables us to predict the x-dependence so that one has 
to rely on parametrizations which are fitted to the data. A more detailed description of the 
structure functions is provided by the parton model which can be applied if one can neglect 
power corrections of the type (1/Q 2 ) P (higher twist effects). Here one asumes that in the Bjorken 
limit (Q 2 — >• 00, x is fixed) the interaction between the hadron and lepton in process (1) proceeds 
via the partons (here the quarks and the gluons) of the hadron. If the scattering of the lepton 
with the partons becomes incoherent the structure function can be written as 
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with similar expressions for the twist two contributions to the spin structure functions g\(x, Q 2 ) 
in which case we introduce the notations AE, AG, AC^i etc.. The vector- and axial-vector 
electroweak couplings of the standard model are given by and afjp respectively with V = 

^y,Z,W and k = l(u), 2(d), 3(s) Further rif denotes the number of light flavours and // 

stands for the factorization/renormalization scale. The singlet (E) and non-singlet combinations 
of parton densities (A fe , Vt) are defined by 



e(^ 2 ) = -E(/*(^VM^ 2 )) 

n / fc=l 

A fc (z,/i 2 ) = f k (z,fi 2 ) + h(z,fi 2 ) - E(z, // 



(7) 
(8) 



V fc (z, /i 2 ) = / fc (z, /i 2 ) - / fi (z, /i 2 ) (9) 

where /fc,/*: denote the quark and anti-quark densities of species k respectively. The gluon 
density is defined by G(z,fi 2 ). The same nomenclature holds for the coefficient functions 
Ci,i(l = q,g) which can also be distinguished in a singlet (S) and a non-singlet (NS) part. 
Like in the case of the structure functions the x-dependence of the parton densities cannot be 
determined by perturbative QCD and it has to be obtained by fitting the parton densities to the 
data. Fortunately these densities are process independent and they are therefore universal. This 
property is not changed after including QCD radiative corrections. It means that the same par- 
ton densities also show up in other so called hard processes like jet production in hadron-hadron 
collisions, direct photon production, heavy flavour production, Drell-Yan process etc. Another 
firm prediction of QCD is that the scale (//) evolution of the parton densities is determined by 
the DGLAP [2] splitting funct ions Pij (i,j — q,q,g) which can be calculated order by order in 
the strong coupling constant a s . The perturbation series of gets the form 

P kl = a s PS ) + a 2 p£> + a 3 s P% ] + . .. (10) 

with a s = a s (/! 2 )/47r. The splitting functions P^ are related to the anomalous dimensions 7^ 
corresponding to twist two local operators Of 1 '"' ln (x) of spin n via the Mellin transform 

7§° = - f dzz"' 1 P^z) (11) 

J o 
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These operators appear in the light cone expansion of the product of two electroweak currents 
which shows up in the calculation of the cross section of process (1) 

oo 

j(x) j(o) -EE of -"-(o, /i 2 ) (12) 

n=0 fe 

where (12) are the Fourier transforms of the coefficient functions cj^ (5), (6) (k = q,g) in 
Minkowski space (x M ). Like the splitting functions they are calculable order by order in a s and 
the perturbation series takes the form 

Q, k = 6 kq + a s Cg + a 2 s eg + a 3 s cfj + ... (13) 

with i = 1,2,3, L and k — q,g. We will now review the higher order QCD corrections to the 
splitting functions and the coefficient functions which have been calculated till now. 



2 Splitting Functions 

The splitting functions are calculated by 

1. Pip Gross and Wilczek (1974) [3]; Altarelli and Parisi (1977) [2]. 

2. AP/f Sasaki (1975) [4]; Ahmed and Ross (1976) [5]; Altarelli and Parisi [2]. 

3. P^p Floratos, Ross, Sachrajda (1977) [6]; Gonzales-Arroyo, Lopez, Yndurain (1979) 

[7]; Floratos, kounnas, Lacaze (1981) [8]; Curci, Furmanski, Petronzio (1980) [9]. 

4. AP$p Zijlstra and van Neerven (1993) [10]; Mertig and van Neerven (1995) [11]; 

Vogelsang (1995) [12]. 

Notice that till 1992 there was a discrepancy for P^ between the covariant gauge [6-8] and 
the lightlike axial gauge calculation [9] which was decided in favour of the latter by Hamberg 
and van Neerven who repeated the covariant gauge calculation in [12]. The DGLAP splitting 
functions satisfy some special relations. The most interesting one is the so called supersymmetric 
relation which holds in M = 1 supersymmetry [13]. Here the colour factors, which in SU(N) 
are given by C F = (N 2 - 1)/2N , C A = N, T f = 1/2 become C F = C A = 27) = N. The 
supersymmetric relation then reads 

pS,(k) , p(k) _ p(k) _ p(k) _ n / 14 \ 
<2<? ~ 91 1 qg 1 gg w V- 1 ^; 



A pS,(k) + A p(k) _ A p(k) _ A p(k) = Q / 15 N 

qg ^ gq qg gg ^ ' 

which is now confirmed up to first (k = 0) and second (k = 1) order in perturbation theory. The 

(2) a (2) 

third order splitting functions P^ ,AP^' are not known yet. However the first few moments 
^2),(n) n _ 2 ? 4 ; g ; 8, 10 have been calculated by Larin, van Ritbergen, Vermaseren (1994) [14]. 
Besides exact calculations one has also determined the splitting functions and the anomalous 
dimensions in some special limits. Examples are the large rif expansion carried out by Gracey 
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(1994) [15]. Here one has computed the coefficients 6 2 i and b 3i in the perturbations series of the 
non-singlet anomalous dimension 



Iqq 



a 2 [n f C F b 2 i] + a 3 s n 2 f C F b 31 



+ n f C A C F b 32 + n f C F b. 



':« 



+ ... 



(16) 



Further Catani and Hautmann (1993) [16] calculated the splitting functions Pij(x) in the limit 
x — > 0. The latter take the following form 



PPM 



In x 



Ak),(n) 
lij 
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(17) 



The above expressions follow from the BFKL equation [17] and /^-factorization [18]. Some 
results are listed below. The leading terms in are given by 
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where £(n) denotes the Riemann zeta-function. Further we have in leading order l/(n — 1) 
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Kirschner and Lipatov (1983) and Blumlein and Vogt (1996) have also determined the subleading 
terms in the splitting functions (anomalous dimensions). They behave like 



ij 



ln 2k z 



lij 



n— 



2k+l 



n 



(22) 



The same logarithmic behaviour also shows up in AP^- and A 'yff- In the latter case the 
expressions in (22) become the leading ones since the most singular terms in (17) decouple in 
the spin quantities. The expressions in (22) have been calculated for the spin case by Bartels, 
Ermolaev, Ryskin (1995) [20] and by Blumlein and Vogt (1996) [21] who also investigated the 
effect of these type of corrections on the spin structure function gi(x, Q 2 ). Finally the three-loop 
anomalous dimension is also known (see Chetyrkin,Kuhn (1993) [22] and Larin (1993) 

[23]). It reads 



A 7 ? (1) 



iqq 



a 2 s [-6n, C F \ + a 3 s 



18CI - 



142^ „ \ 4 2 ' 
— C A C F )n f + - n f C F 



(23) 



Notice that the second order coefficient was already determined by Kodaira (1980) [24]. 
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3 Coefficient Functions 



The higher order corrections to the coefficient functions are calculated by 

1. , CQ i = 1,2,3, L Bardeen, Buras, Muta, Duke (1978) [25], 

see also Altarelli (1980) [26]. 

2. AC^ , AC« Kodairaet al. (1979) [27], 

see also Anselmino, Efremov, Leader (1995) [28] 

Together with the splitting functions Pjj \ AP^ (k — 0, 1) one is now able to make a complete 
next-to-leading (NLO) analysis of the structure functions Q 2 ) (i = 1, 2, 3, L) and g\{x, Q 2 ). 
The second order contributions to the coefficient functions are also known 

1. C$,C?j i = 1, 2, 3, L Zijlstra and van Neerven (1991) [29] 

2. ACj 2 >, ACj 2 > Zijlstra and van Neerven (1993) [10] 

The first few moments of (i = 2, L; k = q, g) were calculated by Larin and Vermaseren 
(1991) [30] and they agree with Zijlstra and van Neerven [29]. The first moment of AC^ was 
checked by Larin (1993) [31] and it agrees with the result of Zijlstra and van Neerven [10]. The 
third order contributions to the coefficient functions are not known except for some few moments. 
They are given by 

1. C{ 3) q {1) (Bjorken sum rule) Larin, Tkachov, Vermaseren (1991) [32]; 

2. (Gross-Llewellyn Smith sum rule) Larin and Vermaseren (1991) [33]; 

3. AC^'^ (Bjorken sum rule ) Larin and Vermaseren (1991) [33]; 

4. Cf g ),(n) (i — 2, L) n = 2, 4, 6, 8 Larin, van Ritbergen, Vermaseren (1994) [14], 

(see also [34]). 

(2) (2) 

Since the three-loop splitting functions P^ , AP^ are not known, except for a few mo- 
ments, it is not possible to obtain a full next-to-next-to-leading order (NNLO) expression for the 
structure functions. However recently Kataev et al. (1996) [35] made a NNLO analysis of the 
structure functions F 2 (x, Q 2 ), F 3 (x,Q 2 ) (neutrino scattering) in the kinematical region x > 0.1 
which is based on 7 ^( 2 )'(") for n = 2,4,6,8, 10 [14]. Like in the case of the DGLAP splitting 
functions Catani and Hautmann (1994) [16] also derived the small x-behaviour of the coefficient 
functions. At small x the latter behave like 
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[3mm] The ingredients of the derivation are again the BFKL equation [17] and /c T -factorizaton 
[18]. from [16] we infer the following Mellin-transformed coefficient functions. 
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The order a 2 coefficients were already obtained via the exact calculation performed by Zijlstra 
and van Neerven (1991) [29]. The subleading terms given by 



2.k 



ln 2 '- 1 ^ 



(/>!) 



(29) 



were investigated by Bliimlein and Vogt (1996) [21]. The most singular terms shown in (24) 
do not appear in the spin coefficient functions AC^ because the Lipatov pomeron decouples in 
polarized lepton-hadron scattering. Therefore the most singular behaviour near x = is given by 
(29) (see [10], [21]). Besides the logarithmical enhanced terms which are characteristic of the low 
x-regime we also find similar type of logarithms near x — 1. Their origin however is completely 
different from the one determining the small a;-behaviour. The logarithmical enhanced terms 
near x — 1, which are actual distributions, originate from soft gluon radiation. They dominate 
the structure functions Fj and ^ near x — 1 because other production mechanisms are completely 
suppressed due to limited phase space. Following the work in [36] and [37] the DGLAP splitting 
functions and the coefficient functions behave near x — 1 like 



pNS,(k) = a pNS,(k) 

gg gg 



1 



x 



p(n) = ^p(k) ^ 
99 99 
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X 



AC NS,(k) 



c 



NS,(k) 



In 



2k-l, 



l-x) 



■t,q 



X 



{1 = 1,2,3) 



(30) 



(31) 



Notice that the above corrections cannot be observed in the kinematical region (x < 0.4) acces- 
sible at HERA. Furthermore the behaviour in (30) is a conjecture (see [7]) which is confirmed 
by the existing calculations carried out up to order a 2 s . 



4 Heavy Quark Coefficient Functions 

The heavy quark coefficient functions have been calculated by 



1. C«(a;,Q 2 ,m 2 ) ( 

2. ACW^Q 2 ,™ 2 ) 



2,L) 



Witten (1976) 



Vogelsang (1991) [39]; 



3. cg(x,Q 2 ,m 2 ), C%(x,Q 2 ,m 2 ) (i = 2,L) 



(2), 



Laenen, Riemersma, Smith, van Neerven 
(1992) [40]. 
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where m denotes the mass of the heavy quark. The second order heavy quark spin coefficient 
functions AC^(x,Q 2 ,m 2 ) and AC^ 2 \x,Q 2 ,m 2 ) are not known yet. Due to the presence of the 
heavy quark mass one was not able to give explicit analytical expressions for C it k(i — 2, L; k — 
q,g). However for experimental and phenomenological use they were presented in the form 
of tables in a computer program [41]. Analytical expressions do exist when either x — > or 
Q 2 3> m 2 . In the former case Catani, Ciafaloni and Hautmann [42] derived the general form 



r (i) 



~ - W- 2 (x) f(Q 2 , m 2 ) (/ > 2,i = 2, L; k = q, g) (32) 

x >u x 



Like for the light parton coefficient functions (see (24)) the above expression is based on the 
BFKL equation [17] and /^-factorization [18]. In second order Buza et al. (1996) [43] were able 
to present analytical formulae for the heavy quark coefficient functions in the asymptotic limit 
Q 2 3> m 2 . This derivation is based on the operator product expansion and mass factorization. 



5 Phenomenology at low x 

Since the calculation of the higher order corrections to the DGLAP splitting functions P^- and 
the coefficient functions Cik is very cumbersome various groups have tried to make an estimate 
of the NNLO corrections to structure functions in particular to F 2 (x, Q 2 ). The most of these 
estimates concerns the small rc-behaviour. In [44] Ellis, Kunszt and Levin Hautman made a 
detailed study of the (Revolution of P 2 using the small x-approximation for P^ (17) and C^l 
(24). Their results heavily depend on the set of parton densities used and the non leading small x- 

(2) 

contributions to . The latter are e.g. needed to satisfy the momentum conservation sum rule 
condition. Large corrections appear when for x — > the gluon density behaves like xG(x,fi 2 ) — > 
const, whereas they are small when the latter has the behaviour xG(x, fi 2 ) — >■ x~ x (\ ~ 0.3 — 0.5 
; Lipatov pomeron). 

However other investigations reveal that the singular terms at x — 0, present in P^ and C^k, do 
not dominate the radiatve corrections to F2(x,Q 2 ) near low x. This became apparent after the 
exact coefficient functions or DGLAP splitting functions were calculated. 

In [45] Gliick, Reya and Stratmann (1994) investigated the singular behaviour of the second 
order heavy quark coefficient functions (32) in electroproduction and they found that its effect 
on F 2 was small. 

Similar work was done by Blumlein and Vogt (1996) [21] on the effect of the logarithmical terms 
(22), (29) on gi(x,Q 2 ) which contribution to the latter turned out to be negligable. 
Finally we would like to illustrate the effect of the small x-terms, appearing in the coefficient 
functions C^l and C^ k , on the structure functions F 2 (x, Q 2 ) and Fl(x, Q 2 ). For that purpose we 
compute the order a 2 contributions to F 2 and Fl. Let us introduce the following notations. When 
the exact expressions for the coefficient functions are adopted the order a 2 contributions to 
Pj will be called §F^' exact . if we replace the exact coefficient functions by their most singular 
part which is proportional to \jx (see (24)) the order a 2 contributions to Pj are denoted by 
8F^' app . The results are listed in table 1 and 2 below. Further we have used the parton density 



7 



sets MRS(DO) {xG{x,fi 2 ) -»• const, for x -> 0) and MRS(D-) (xG(x, u 2 ) -»• ar A for x -> 0) [46] 







MRS(DO) 






MRS(D-) 






| F 2 NLO 


<• t-i(2), exact 
Of 2 


5F 2 {2) ' app 


F 2 NLO 


r t-i(2), exact 
Of 2 


5F 2 (2) ' app 


io- 3 


0.67 


-0.069 


0.088 


0.99 


-0.084 


0.116 


10" 4 


0.82 


-0.088 


0.158 


2.29 


-0.226 


0.349 


10~ 5 


1.00 


-0.092 


0.251 


5.99 


-0.665 


1.059 






^p(2),exact 
L 






^ p{2) , exact 


fip<j?)< a PP 


10" 3 


0.149 


-0.029 


-0.040 


0.263 


0.008 


-0.052 


10" 4 


0.210 


-0.062 


-0.071 


0.780 


0.031 


-0.156 


10~ 5 


0.281 


-0.102 


-0.113 


2.370 


0.105 


-0.475 



From the table above we infer that a steeply rising gluon density near x = (MRS(D-)) 
leads to small corrections to F 2 and F L . On the other hand if one has a flat gluon density 
(MRS(DO)) the corrections are much larger in particular for F L . A similar observation was 
made for F2 in [44]. However the most important observation is that the most singular part of 
the coefficient functions gives the wrong prediction for the order a 2 contributions to the structure 
functions except for Fl provided the set MRS (DO) is chosen. This means that the subleading 
terms are important and they cannot be neglected. Therefore our main conclusion is that only 
exact calculations provide us with the correct NNLO analysis of the structure functions. The 
asymptotic expressions obtained in the limits x — > 0, x — > 1 and Q 2 3> m 2 can only serve as a 
check on the exact calculations of the DGLAP splitting functions and the coefficient functions. 
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